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Abstract 

Stress and strain are important factors in determining the mechanical, electronic, and optical 
properties of materials, relating to each other by the material's elasticity or stiffness. Both are 
represented by second rank field tensors with, in general, six independent components [T]. Mea- 
surements of these quantities are usually achieved by measuring a property that depends on the 
translational symmetry and periodicity of the crystal lattice, such as optical phonon energies using 
Raman spectroscopy P], the electronic band gap using cathodoluminescence [3], photoelasticity 
via the optical birefringence 0], or Electron Back Scattering Diffraction (EBSD) [5]. A recipro- 
cal relationship therefore exists between the maximum sensitivity of the measurements and the 
spatial resolution. Furthermore, of these techniques, only EBSD and off-axis Raman spectroscopy 
allow measurement of all six components of the stress tensor, but neither is able to provide full 
3D maps. Here we demonstrate a method for measuring the full stress tensor in diamond, using 
the spectral and optical polarization properties of the photoluminescence from individual nitrogen 
vacancy (NV) colour centres. We demonstrate a sensitivity of order 10 MPa, limited by local 
fluctuations in the stress in the sample, and corresponding to a strain of about 10~^, comparable 
with the best sensitivity provided by other techniques. By using the colour centres as built-in 
local sensors, the technique overcomes the reciprocal relationship between spatial resolution and 
sensitivity and offers the potential for measuring strains as small as 10~^ at spatial resolution of 
order 10 nm. Furthermore it provides a straightforward route to volumetric stress mapping. Aside 
from its value in understanding strain distributions in diamond, this new approach to stress and 
strain measurement could be adapted for use in micro or nanoscale sensors. 



The NV centre is the most common colour centre in diamond, comprising a substitutional nitrogen 
atom and a vacancy in an adjacent site along one of the four (111) crystal axes, as depicted in figure 
n[a) . The perturbation of its optical emission by applied uniaxial stress via the piezospectroscopic effect 
[3 was characterised in the seminal work of Davies and Hamer in 1976 [8^ and used to identify the 
centre as possessing C^y symmetry and a. E ^ A optical transition. In recent years, optical microscopy 
of individual centres has become popular since they were demonstrated to be exceptionally stable single 
photon emitters [S] and to possess an electronic structure that couples the highly coherent spin state of 
the centre to its optical transitions, |10j thus offering considerable potential for quantum information 



1 



applications. As a result of this renewed interest, the electronic structure of the colour centre has 
been studied extensively and is now well understood in terms of the orbital energies, spin-orbit and 
spin-spin interactions [111 1121 [T3] . and hyperfine coupling. The perturbation of the energy levels by 
external electric and magnetic fields is also established, leading to demonstrations of the use of NV 
centres as nanoscale sensors [HI I28j . 




Figure 1: The NV centre structure: (a) Schematic of the NV structure and definition of NV axes 
(x, y, z) and crystal axes {X, Y, Z). The cube represents | of the cubic unit cell of the diamond crystal. 
The nitrogen atom is colored green and the vacancy is represented by the dashed circle in the center 
of the cube. The three carbon atoms adjoining the vacancy are in grey, (b) Rotation of the E^' and 
Eyi orbitals relative to NV axes x and y under a symmetry-breaking stress. E^' is taken to be the 
higher energy orbital, consistent with (c) Projection onto the (001) image plane (XY) of the four 
possible orientations of NVs relative to the crystal axes, (d) Projection of the E^'^y' orbitals in (b) 
onto the image plane, leading to the measured data in figure [3] 

The effect of local stress on a single NV centre can be understood by considering the effect that 
different 'components' of the stress field have on the defect symmetry. To this end we define a set of 
orthogonal axes (x, y, z) with z parallel to the NV axis and x lying in one of the three mirror planes. 
For example an NV centre with its (z) axis along [111] we choose x to lie along [112] and y along [110] 
as shown in figure [ija). 

The ^E optically excited state of the NV centre has twofold spatial degeneracy and threefold spin 
degeneracy, and under the C^y symmetry of the defect these six eigenstates are split into four levels by 
the spin-orbit and spin-spin interactions, with splittings of order a few GHz. Under stress the energies 
of the ^E states are perturbed. A stress that preserves the C31, symmetry will result in a global shift 
of the ^E levels, whilst a stress that breaks the symmetry modifies the splittings in the ^E manifold 
and rotates the orbitals in the (x, y) plane. In this respect the stress perturbation is similar to that of 
a vector with three components {£,XT^y,£,z) [IHlllH]- In the limit that the stress perturbation greatly 
exceeds the spin-orbit and spin-spin splittings the six '^E states form two branches corresponding to 
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orthogonal spatial orbitals E^' and Eyi , and the perturbation can be expressed as a secular matrix in 
the {Ex,Ey) basis [T^ . 



+ /3 7 
7 a — P 



(1) 



where a, /3 and 7 are energy shifts determined by the six components of the stress tensor. The 
eigenvalues of this secular matrix are the energy changes of the stress-shifted E states, and the eigen- 
vectors provide the orientations of x' and y', the '^E -^^ A transition dipoles, in the {x,y) plane. The 
optical transition line therefore becomes a doublet with a mean energy shift of a and a splitting of 
S — 2-y/ -f 72^ and the E dipoles are rotated by an angle (p = 1/2 arctan(7//?) relative to x and y 
(figure [if^b)). For a [111] oriented NV centre the energy shifts are given by 



a = Ai (cr^^ + a^^ +<yzz) + 2^2 {(^yz + ^ zx + ^ xy) 

I3 = B {2a - cr^^ - cr^^ ) + C (2cr^^ - cr^^ - a^^ ) (2) 
7 = (ct^^ - cr^,.) + \/3C (0-^2 - CT^x) 

where the parameters Ai, A2, B, and C have been determined as 1.47, -3.85, -1.04, and -1.69 
meV/GPa respectively [5], and the stress components aij (where i,j = X,Y,Z) are referred to the 
crystal axes. The equations for an NV centre with one of the three other orientations are obtained by 
performing simple rotations on the coordinate axes (X, Y, Z) (see supplementary information) . The 
equivalent equations to ^ with the stress referred to the NV axes are 

a = ^1 (^.. + ^yy + '^.J + A2 (2(7„ - a_ - a J 

/3 = iB + C){a^^-aJ+V2i2B-C)a^^ (3) 
7 = -2iB + C)a^^ + V2{2B - C)a^^^ 

from which it is clear that the global shift a results from any stress that shares the NV symmetry 
(^z), /3 results from stress that shares only the xz mirror plane symmetry (S,x)^ and 7 results from 
stress with neither the NVs cylindrical nor its xz mirror symmetry (^y). 

Since the stress tensor contains six independent parameters, measurement of a single NV centre 
can clearly provide only partial information of the local stress field. To obtain all six components of 
the tensor it is therefore necessary to measure three NV centres that experience the same stress field 
but which have different orientations in the lattice. With a, /3 and 7 determined for three differently 
oriented colour centres one obtains an overdetermined set of nine linear equations with six unknown 
variables, from which a unique solution is obtained using the Moore- Penrose pseudoinverse method 

To demonstrate this method we present photoluminescence data from isolated NV centres in a 
single crystal diamond grown by plasma assisted chemical vapour deposition. A sample was selected 



which displays isolated NV centres in a region of particularly high grown- in stress (figure 2(a) ), so that 
at a temperature of 77 K, many of the zero phonon lines reveal clear splittings of up to 0.8 nm or 2.4 
meV (figure [2(b)[ ). The viewing axis is Z || [001], such that the integrated emission of a given NV 
centre is approximately 3 : 1 polarized, allowing easy discrimination between NVs with axis parallel 
to [111] or [111] and those with axis parallel to [111] or [111] (figure [l]^c)). The ambiguity (effectively 
a mirror in Z) in identifying the orientation of a given NV centre is not possible to resolve with this 
sample geometry, but as we shall see presents only a minor limitation to the results. 

The two E ^ A electronic transitions emit linearly polarized photons aligned with the dipole axes. 
The measured projections of these polarizations onto the (001) image plane allow determination of the 
angle (/> and thus we obtain 7 = 5sin(2(/)) and 13 = S cos(2(/)). The parameter a is taken as the energy 
shift between the mid point of the two lines and 1.945 eV (637.5 nm). 

As the z-projections of individual NVs are not known we performed measurements on several 
proximal centres and looked for groupings of parameters to identify different defect orientations. An 
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defect 


angle 


a{meV) 


l3{meV) 


j{meV) 


A 


90 


1.966 


1.712 


-0.688 


B 


90 


2.359 


1.391 


-0.563 


C 


90 


2.354 


1.347 


-0.599 


D 





1.322 


0.495 


1.376 


E 





1.137 


0.163 


1.191 


F 





1.694 


0.304 


1.087 


G 





1.740 


0.354 


1.052 



Table 1: Polarisation and energy shift data for seven NV centres labeled in figure |2j 



example of a group of seven NVs studied is shown in figure 2(a) Polar plots revealing the optical 
polarization from each of these centres are shown in figure [3j and the corresponding a, /3, 7 values 
listed in table [T] 




Figure 2: Low temperature Photoluminescence from individual nitrogen vacancy centres in highly 
strained material: (a) PL image of the sample area of interest with the seven NVs used in this study 
labeled A-G; (b) PL spectrum of a single NV centre with a split zero phonon line (inset). 



It can be seen from the tabulated data that in two cases pairs of NVs that are spatially close (B & 
C, and F & G) give very similar values for all parameters measured. The conclusion is that these NVs 
are co-oriented, and are not affected substantially by local perturbations to the lattice. Of the 22 NVs 
measured in total for this study, we found six examples of pairs of proximal NVs whose parameters 
agreed to better than 10%, suggesting that for these NVs the assumption of a dominant long-range 
stress field is valid. Several NVs also revealed behaviour that was inconsistent with the assumption, 
presumably resulting from the presence of a strong local stress field due to a nearby defect such as a 
lattice vacancy, or an electric field due to a trapped charge. 

We now use the data for NVs A, (B or C) and D to calculate the local stress tensor in the top right 
of the image. Choosing the orientations A|| [111], (B or C)|| [111] and D|| [111] yields the stress tensor 

290(5) 90(0) ±6(0) 

90(0) 758(0) ±45(0) | MPa (4) 
±6(0) ±45(0) -14(2) 
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Figure 3: Measured polarisation data for NVs A, B, C and D. The radius represents the emission 
intensity and the angle is that of the polarizing filter, The solid curves are fits to the function 
Iq + AIcos'^{(j)p + (p') where Iq, AI, and (p' are fitting parameters. 



where the figures presented are the mean values (standard deviations) for the tensor components 
calculated using the sets {A,B,D} and {A,C,D}. The ± values in axz and ayz indicate the effect 
of inverting Z for all three centres. Reversing the choice of orientation for D (only) creates a distinct 
configuration and results in the tensor 

259(1) 98(2) ±44(2) \ 

98(2) 758(4) ±27(2) MPa (5) 
±44(2) ±27(2) -31(6) / 

The most significant differences between a and cr' are in the tensor components containing Z, 
corresponding to the Z reflection of centre D. We note however that the choice of NV B or C makes very 
little difference to the results in either case, demonstrating that the difference in the stress experienced 
between these two NVs is of the order of a few MPa at most. That the standard deviations are 
consistently lower in a than in a' may further suggest that this choice of orientation for D is the 
correct one of the two. 

From a calculated stress tensor we can use equations ^ to recalculate a, (3 and 7 for each of the 
NV orientations thus providing a convenient test of the consistency of the measurement. Tensor a 
yields the results presented in table [2j 

From these results we can clearly identify the [111] values as relating to NV D, [111] relating to NVs 
B and C, and [111] relating to NV A. The discrepancies between the recalculated values and measured 
values are of order 10%, consistent with the observed variation in parameters described previously. 
Noteworthy is the absence of any NVs with parameters comparable to the predictions for the [ill] 
orientation. We attribute the orientations of all three (more distant) NVs E, F and G as [111] the 
differences between their parameters and those of NV D being due to longer range changes in the stress 
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NV axis 


angle 


a{m,eV) 


/3{meV) 


j{meV) 


[111] 





1.219 


0.729 


0.957 


[111] 


90 


2.511 


1.358 


-0.993 


[111] 





0.434 


0.901 


0.731 


[111] 


90 


1.914 


1.489 


-0.695 



Table 2: Recalculated energy perturbations for the four NV orientations using stress tensor a. 
field. 

Discussion 

The tensors in equations Q and ([5| indicate a stress that is primarily biaxial in the (001) plane, with 
principal axes closely aligned to the [100] and [010] crystal axes. That there is very little tensile or 
compressive stress along the [001] axis is not surprising as this is the growth axis of the crystal and the 
NV defects measured here are ~ 50 |^m below the surface. The apparent absence of any NVs oriented 
parallel to [111] may be a result of both the stressed lattice and the NV defects having been generated 
during the CVD grown process, whereupon there are likely to be preferred and inhibited directions for 
defect orientation }21j . 

We now discuss issues related to the sensitivity of the technique. For simplicity we have presented 
only the theory applicable to a stress perturbation larger than 10 MPa, equivalent to the ~10 GHz 
splitting of the unperturbed '^E level, and have treated the spin levels as degenerate. Smaller stresses 
can equally be determined by measuring in more detail the spectrum of the state and using the full 
theoretical description that includes the spin-orbit and spin-spin interactions jl2l I13j . The sensitivity is 
also clearly limited by the spectral resolution of the measurements, which is ultimately determined by 
the transition line widths. These are a strong function of temperature, so that the sensitivity increases 
as the sample is cooled. Line widths approaching 10 MHz have been observed using photoluminescence 
excitation spectroscopy at 2.8 K [32], corresponding to stress magnitudes of order kPa, or strains of 
order 10~^. Although at such extreme sensitivity local variations in the built-in lattice stress would 
certainly dominate, changes in strain of this magnitude may ultimately be detectable. At higher 
temperatures the line widths broaden rapidly due to the dynamic Jahn Teller effect [23] , and at 77 K 
they are limited to about 10 GHz corresponding to stresses of order 10 MPa. At 300 K the zero phonon 
line is generally not clearly visible, and so stress mapping under ambient conditions would require 
utilisation of the stress perturbation to the electron spin state of the centre, measured using optically 
detected magnetic resonance (ODMR). The spin-orbit interaction introduces a stress component to the 
spin splittings in both the orbital doublet excited state 'ST and, through inter level mixing, the orbital 
singlet ground state |28) . Although the stress dependence of the spin transitions is much stronger in the 
excited state, the ODMR line width is limited by the optical relaxation time (~12 ns) thereby limiting 
the sensitivity of excited state techniques. Ground state ODMR probing the spin transition at 2.88 
GHz is therefore likely to be the more sensitive method. In this way, Dolde et al have demonstrated 
the use of single NV centres as electric field sensors with a dc sensitivity of 891 ±20 V cm^^ Hz~2^ 
equivalent to a stress of about 90 kPa Hz~ s . By combining data from three differently oriented NV 
centes the full stress tensor can similarly be extracted (see supplementary info for details). 

Volumetric imaging of the stress tensor is possible by recording a 3-D distribution of defects. With 
a standard scanning confocal optical microscope this can provide sub-micrometre spatial resolution 
per defect, therefore up to ~ 1 |i.m resolution for the full tensor. However NV centres in diamond 
are also ideally suited to super-resolution optical techniques such as stimulated emission depletion 
microscopy fB] that utilise an intense depopulating beam and require highly robust photon emitters. 
Such techniques are entirely compatible with the methodology described here and could yield spatial 
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resolution on a length scale of a few nanometres for systematic investigations of local stress fields. 
We note that such spatial resolution would be fundamentally incompatible with the various existing 
approaches to stress/strain measurement. EBSD, Raman, and band gap luminescence all interrogate 
quantities that depend on the periodicity of the crystal lattice and can therefore only measure strain 
averaged over a volume, or length, that is large compared with the lattice spacing. The use of point 
defects as sensors therefore opens up new possibilities for stress/strain measurement such as measuring 
lattice scale fluctuations in stress caused by nearby defects or surfaces, or for use in nanoscale sensors, 
for example in nano-electromechanical systems (NEMS). 



The sample studied was a high purity, [001] oriented type Ila diamond grown by microwave plasma- 
assisted chemical vapour deposition (CVD) [T7]. Homoepitaxial CVD diamond growth was performed 
using a selected high pressure high temperature grown substrate, which had been carefully processed 
prior to growth in order to minimize the dislocation density in the grown CVD diamond. The resultant 
diamond had a nitrogen concentration of < 5 ppb as measured by ensemble electron paramagnetic reso- 
nance (EPR) and a boron concentration of < 0.5 ppb as measured by secondary ion mass spectrometry 
(SIMS). After growth the sample was mechanically polished to have a surface RMS roughness < 1 nm 
over a 1 X 1 |j.m area. 

Low temperature photoluminescence microscopy and spectroscopy was recorded using a custom- 
made fiber-coupled scanning confocal microscope |18j . Photon counting for image generation utilised 
a Perkin Elmer SPCM. Spectra were recorded using an Acton SP 500 monochromator and Princeton 
Spec 10:100B CCD camera. Polarization dependence of emission was measured by placing a linear 
optical polarizer in front of the collection aperture of the microscope. A slight wedge in the plate leads 
to partial misalignment of the microscope as the polarizer is rotated, resulting in asymmetry in some 
of the data in figure |3] 

A Supplementary information 

A.l a, (3 and 7 for all NV orientations 

The expressions for a, (3 and 7 depend on the orientation of the NV centre. Since the NV axis can 
take one of the four < 111 > directions (the order of the nitrogen and vacancy sites is not important) 
we require four sets of equations. To calculate these we start with the equations for the [111] oriented 
defect with the stress tensor components referred to the crystal axes A", Y ^ and Z (see equation ([2])) 
and apply successive rotations of ^ around the z axis to the coefficients matrix M = {cTij}. Since 
the stress component Uij transforms under symmetry operations as the product ij , the result of these 
rotations can be seen by the successive mappings X ^ Y and Y 1— > —AT. The full set of equations is 
presented below. 

For the [111] oriented NV centre: 



Methods 



13 = B {2a - a^^ -a^^) + C {2a - a^^ - cr^x ) 
7 = \/3B (cr^^ - cr^^) + \/3C(o-^2 -<^zx)- 



(6) 



For a [111] oriented NV centre: 



a = Ay {a^^ + a^^ + cr^z) + 2^2 (-o-^z + 
/3 = B {2a - a^^ -a^^) + C {-2a + a^^ 
7 = Vl5B {-a^^ + cr^^y- ) + {a^^ + cr^x ) • 




(7) 
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For the [111] oriented NV centre: 



P^B {2a,, - -G^^) + C {2g^^ + G^, + G,^ ) (8) 

7 = ^/?>B (cr^^ ~(7yy) + \/3C (-0-^2 + C^zx ) ■ 

For a [111] oriented NV centre: 

Ai (o-^^ + (Tyy + (Tzz) + 2yl2 (+0-^2 - g,^ - g^^) 
(i = B {2g,, - G^^ -G^^) + C {-2g^^ - G^, + G,^ ) (9) 



A. 2 The effect of different viewing orientations 

The use of a [001] oriented sample in this study derives from the standard growth conditions of high 
purity CVD diamond and the ready availabihty of a suitable sample for study. However as mentioned 
earlier this viewing axis results in an unavoidable ambiguity in the orientation of any given NV centre, 
effectively a mirror in Z, so that for a given NV centre we can not distinguish between [111] and [111], 
or between [Til] and [111] orientations. As discussed in the text, we can in some cases make good 
estimates of which NVs share the same orientation based on their common behaviour in a uniform 
stress field, but it is not possible to identify unambiguously what that orientation is. There are two 
ways of avoiding this problem and thereby providing unambiguous identification of the NV orientations 
and thus the full stress tensor, subject only to the uniformity of the stress field. The first is to choose 
a viewing axis that provides straightforward identification of the NV orientation. For example, a [111] 



viewing axis results in the projections of the orbital planes onto the focal plane as shown in figure 4(a) 
with four different 'ellipses' for the four different NV orientations. The result is that the dependence 
of PL intensity on excitation polarisation allows simple determination of the orientation of each NV 
centre. 

[Ill] 



[111] 



[111] 




[111] [111] 



(a) view along [111] direction (b) view along [110] direction 

Figure 4: Projections of the orbital planes of the four different NV orientations when viewed along (a) 
[111] and (b) [111] direction respectively. 

Conversely, other viewing axes are less conducive to the stress tensor evaluation. As shown in figure 
|4(b)[ a [110] orientation results in a twofold ambiguity in orientation as does [001]. But, more seriously, 
it does not allow for measurement of the angle of rotation of the E^i , Eyi eigenstates away from the 
x, y axes for two of the orientations since the orbital plane is viewed edge-on. A [110] viewing axis can 
not therefore be used to determine the full stress tensor by this method. 
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The second approach is to resolve the ambiguity with a [001] viewing axis, by applying a field which 
breaks the mirror symmetry in Z (or in both [110] and [110]). This could be achieved by the temporary 
application of an electric field, under which the different orientations would be expected to respond 
differently. For instance an electric field aligned along a (111) axis would allow identification of NV 
centres aligned parallel with that axis, as they would display a global shift of the optical transition 
energies but no change in splitting or polarization. 



A. 3 Stress tensor measurement using spin resonance 

The Ai electronic ground state of the NV centre has a non-zero permanent dipole as a result of a small 
amount of mixing with the E symmetry excited states, which gives rise to a stress and electric field 
dependence of the spin transitions. The components of this dipole aligned parallel and perpendicular 
to the NV axis have been measured as Sgajh = 0.35±0.02Hz cm V~^ and d^^/h = 17±3Hz cm V~^ 
|27) . The Hamiltonian for the ground state in the presence of local stress, plus electric and magnetic 
fields m nam] is given by 



1 



Bgs - {hDgs + n.)[^^ - -S{S + 1)] + Ai^.geS x B 



3 ^ " (10) 

- d^gMAsl - sD + n,(5,5, + SyS.,)\ 

where as defined in figure^ z is along the NV axis and x is in one of the three mirror planes. Dgs 
is the spin-spin interaction strength, /z^ is the Bohr magneton, g^. is the electron g-factor. ^x^^y^ and 
IIj are the components of a vector that sums the electric field and the effect of stress on a single NV 
centre, and S^^ Sy and are the triplet spin operators, shown below in the (x, y, z) basis 






S'^ = ;i OOl| Sy^h\ I S';,=:?i(100| (11) 
\0 1 0^ 

which result in spin eigenstates 

|y) = I H k) = I I (12) 





while S{S + 1) = + Sy + = 2. Matrix elements can thus be calculated using equations ( 10 1, (111, 



and (12) to write the Hamiltonian in the (x,y,z) basis 



hDgs + dlsU^ + dj,U^ -dj.Uy 

Hgs = I -dj^^iiy hDgs + 4«n, - 4,n, I (13) 

~2{hDgs+4s^, 

where we have taken the case of zero magnetic field for simplicity (we note that, as pointed out by 
Dolde et al in [28], care must be taken to minimise the magnetic field, especially its z component, to 
achieve high sensitivity in these measurements). The effect of the combined stress and electric field 
on the transitions between the z spin state and the other two eigenstates is therefore described by the 
secular matrix in the spin (x, y) basis 



*^gs + I3gs Igs 
"igs '^gs ~ Pg 



(14) 
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where 

Pgs = 4,n, (15) 

Igs = -dj,X\y 

Equation ( [l4| above can be viewed in direct analogy with the secular matrix for the excited state 
as described in equation ffl. In the absence of an electric field the same symmetry arguments that 
lead to equations ^ and expressing a, fi and 7 for the electric dipole transitions in terms of the 
stress tensor components and coefficients Ai,A2,B and C, also apply here, ags, j3gs and 7gs are there- 
fore related in the same way to the six independent tensor components via four equivalent constants 
^i_gs, ^2_gs, Bgs, and Cgs- This allows us to separate the stress contribution from the electric field 
contribution and write 

ags = Ai_gs (o-^^ + cr„„ + o-^,) + A2_gs {2a, ^ - a^^ - a^J + SdlsE^ 

PgS = {Bgs + Cgs) {<J^^ - <J J + V2{2Bgs " C g s) CJ ^ ^ + d^, iJ, (16) 
= -2{Bgs + Cgs)a^^ + ^2{2Bgs - Cg,)fT^, - dj^Ey. 

For the [111] oriented centre, with stress and electric field now referred to the crystal axes, we have a 
set of three equations analogous to those of equation (2), 

ags = Ai_gs ((Txx + O'l-r + ^zz) + 2A2_gs {cTyz + ^ zx + ^ xy) + ^d^s {E^ + E^ + E^) 

figs = Bgs {2a - a^^ - a^^ ) + Cgs {2a - a^^ ~ <Jzx) + % {'E^ - E^ + 2E, ) (17) 

-igs = V^Bgs {a^^ - CT^l,) + ViCgs (d^^ ' ^ Z x) ' X ~ ^ ^) . 

with the permutations for the other NV orientations as in equation Q. For a single NV centre, a^s, Pgs 
and 7gs can be measured from optically detected magnetic resonance data by the shift and splitting 
in the spin transitions plus the orientation of the transition-inducing microwave field. The same 
pseudoinverse method can thus be applied to determine the full stress tensor from room temperature 
spin resonance data, once the constants Ai gs, A2_gs, Bgg, and Cgs have been established. 
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